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It is known that the interactions of single Skyrmions are asymptotically 
described by a Yukawa dipole potential. Less is known about the interactions of 
solutions of the Skyrme model with higher baryon number. In this paper, it is 
(*h shown that Yukawa multipole theory can be more generally applied to Skyrmion 

interactions, and in particular to the long-range dominant interactions of the 
Q B — 4 solution of the Skyrme model, which models the a-particle. A method 

that gives the quadrupole nature of the interaction a more intuitive meaning in 
_G the pion field colour picture is demonstrated. Numerical methods are employed 

to find the precise strength of quadrupole and octupole interactions. The results 
I are applied to the B = 8 and B = 12 solutions and to the Skyrme crystal. 
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1 Introduction 



50 years ago, T.H.R. Skyrme suggested a nonlinear scalar field theory as a model 
for nucleons and nuclear interactions pQ. While quantum chromodynamics (QCD) 
is currently accepted as the fundamental theory of nuclear interactions, it is nearly 
impossible to work with to model actual nuclei. Therefore, except for very few results 
derived from first principles (e.g. in lattice field theory), to this date, phenomenolog- 
ical theories are mainly used in nuclear physics. The Skyrme model is a field theory 
intermediate betweeen QCD and the point particle quantum mechanics of nucleons. 

After the Skyrme model was first applied by Adkins, Nappi and Witten [2] to 
derive the properties of single nucleons and delta resonances, it has been used to 
derive the properties of several nuclei. The Skyrme model with massive pions has in 
particular been applied to many nuclei with a baryon number divisible by four, in 
which the nuclei are thought to be composed of a number of interacting a-particles 
[3] . The classical B — 4 Skyrmion solution corresponding to the a-particle has a cubic 
shape and is strongly bound. In the case of B = 8, 12 and 32, the B = 4 subunits 
have been found to be in face-to- face contact, to which the methods in this paper can 
be directly applied to find interaction energies. 

1.1 The Skyrme model 

The Skyrme model [I] is a three-dimensional non-linear sigma model. It includes an 
extra term which is necessary for having solitons according to Derrick's theorem. It 
is easiest to write by taking the Skyrme field U(x) to be an SU(2)-valued field. The 
action then takes the form 

S= T2W d^(-^(L M ^) + ^Tr([L ;i ,L 1 ,][^,^])+m 2 Tr(C/-l)) (1) 

where L M = 11^(8^11) and the metric signature used throughout is (+, — , — , — ). This 
is the rescaled form of the theory: classical energy units are 127r 2 F w /(4e) and length 
units 2/ieF^) where is the pion decay constant and e a parameter of the Skyrme 
model. The constant to is the mass of the pion in length units 2/(eF v ). The m-term 
forces U — > 1 as x goes to spatial infinity; in the massless case (to = 0) the same 
is used as a boundary condition. Note that the action in equation ([!]) can also be 
written in terms of the right-current i? M = {d^U)U^ (it has the same form replacing 
L with R). The field equation can be derived [4] to be 

Tr r k (d^ \ L» + ^[L V ,[L V ,L>*]] \ + TO 2 t/) = 0, k = 1,2,3, (2) 

which is also valid for the right current R^. Here Tk are the Pauli matrices. In 
addition to the spatial degrees of symmetry, it is also invariant under the isospin 
rotation U(x) — > AU(x)A~ 1 where A is an SU(2) matrix. 

This model has been suggested by Skyrme [I] as a theory of pions and baryons 
(here: the proton, neutron and delta resonances which can be made from just up and 
down quarks). It is now believed to be a low-energy effective theory for quantum 
chromodynamics, although a formal derivation has not yet been achieved. 

SU(2) is topologically the same as S 3 , as can be seen by writing U = a + in ■ t 
with the restriction a 2 + tt 2 = 1, so that (a, tt) are coordinates on S 3 . Pions are found 
in this model as the perturbative waves of the form U « 1 + iir ■ r, and their mass is 
given by to. Baryons occur as topological solitons. Considering U at any instant of 
time, it is a map U : M. 3 — s- SU(2) = S 3 . As U(oo) = 1 has a definite value, R 3 can be 
compactified using the one-point compactification M 3 U {00} = S 3 . The field U then 
gives rise to a map S 3 — > S 3 that has an associated degree B E Z. In fact there is 
also a (topologically) conserved current 

B» = --^e^ T Tr(L u L a L T ) (3) 
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to which the associated charge, B = J B , is taken as the baryon number of U (and 
is identical to the degree discussed above). Both points of view show that doB = 0. 

B classifies the field configurations U into topological sectors with different baryon 
numbers. Minimising the energy in one sector will achieve a static solution with a 
given baryon number. Quantising the rotational degrees of freedom gives spinning 
solutions which are the true nucleons and nuclei in this model. 



1.2 The B = 1 Skyrmion 

The basic static non- vacuum solution of equation ^ is the so-called hedgehog solution 
[2j|4], which is spherically symmetric in the sense that U(R(A)x) = AUfflA -1 , where 
R(A) is the SO (3) matrix corresponding to the SU(2) matrix A under the canonical 
SU(2) — > SO(3) (two-to-one) homomorphism. The hedgehog solution is believed to 
be the Skyrmion with baryon number one, i.e. the energy minimiser in the B = 1 
sector. The hedgehog solution is of the form 

U{x) =exp(i/(|f|)£-f) (4) 

where x — x/\x\ and f(r) is the profile function. W.l.o.g. /(oo) = so that U(oo) = 1 
(any rotation in real or isospin space of this U would still be spherically symmetric, 
although it would not technically be a hedgehog, which was chosen as a name as all 
the pion fields are pointing outwards). 

By using the spherical symmetry, one can find the differential equation / has to 
obey in order to be a static solution. This is 

(r 2 + 2 sin 2 /) /" + 2rf + sin 2/ (V 2 - 1 - - mV sin / = (5) 

with the boundary conditions /(0) = w and /(oo) = 0. 

Eq. jg) can be solved numerically and the spinning Skyrmions have been analysed 
by Adkins, Nappi and Witten (massless pion case) [2] and Adkins and Nappi (massive 
pions) [5 to find the properties of nucleons and delta resonances. 

Linearising the hedgehog equation ^ yields 

r 2 f" + 2rf - 2/ - m 2 r 2 / = (6) 

which is the modified spherical Bessel equation of first order. This means, considering 
the asymptotic pion fields as massive scalar fields, the hedgehog assumption forces 
an asymptotic dipole field (the same is true in the massless case), given by a rescaled 
first order modified spherical Bessel function of the second kind (which approaches 
when r — > oo) 

f(r) ~ Cke-» + 7) • (7) 

In the massless case, it can actually be shown that for any static Skyrmion the lowest- 
order nonzero multipole is at least a dipole. It cannot have an asymptotic monopole 
[6]. It is not known whether this result holds in the massive pion case, however a 
monopole is ruled out due to symmetry for the hedgehog. 

Actually, there are three massive pion fields. The hedgehog has a dipole falloff in 
each of these components, the dipole moments being 

4ir ( j , 4ir ( Ci j , 4tt (0 j , (8) 

respectively, for the tti, 7T2 and 113 components. They can be put together in one 
matrix, which is a multiple of the identity matrix: 4nCil3. This is really useful 
in case of rotated hedgehogs: Rotated by the SO (3) rotation matrix R, the dipole 
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moment matrix becomes AttC\R. The general asymptotic interaction energy of two 
rotated hedgehogs, 



9tt \\\X\ \X\ 2 \X\ 3 J V / 

+ ^^^ (1) ^ (2) .J , (9) 

is given in [TJ. Here |X| is the distance between the two Skyrmions, R^ and i?' 2 - 1 are 
the rotation matrices (relative to the hedgehog) , and X is the unit vector pointing in 
the direction of separation. 

Using numerics, the constant G\ is easily accessible. The method employed is the 
same as is used later to find asymptotic coefficients in more difficult cases, but in this 
case we only need to solve an ordinary differential equation (ODE). A shooting method 
is applied, choosing a value for C\ (see cq. Q) first. Then, starting from some finite 
point L (L is chosen large so that the asymptotic is valid at L), ODE solution methods 
are used to find the function f(r) numerically, by solving the equation from right to 
leftPj The value at /(0) is then compared with the target value w, and C\ adjusted 
until it is attained with a given precision. Shooting from right to left circumvents the 
difficulty of deciding whether a solution goes to zero at infinity. This is particularly 
difficult for this equation, because missing zero by only a tiny bit, the solution will 
always escape to infinity as r — > oo. 

In figure [T] the value of C\ is plotted against the pion mass parameter m. It is 
quite surprising that C\ has a minimum at positive m. The value C\ — 2.16 at m = 
was known previously [4 . 

2 The interaction energy of Yukawa multipoles 

In this section we derive the expression for the interaction energy of a 2'-pole with 
a 2 m pole in a massive scalar (Yukawa) field theory, needed in the following sections 

1 Using the odeint method from Python's SciPy package [8] 
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Figure 1: Dependence of the hedgehog dipole parameter C± on the pion mass m. 
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to compute Skyrmion interaction energies. Although the expression is simple, it does 
not seem to be well known in this generality. 

A massive scalar field theory with source p is given by the energy 

E = f (UMWit) + \ m2c f + pA d " x ( 10 ) 



" ' 2 
which leads to the linear field equation 

(v 2 - to 2 ) (j) = p . (11) 

The charge distribution of two Yukawa multipoles can be written as 

p(x) = (-l) l Q {1) tl ... H (d n . . . diM*) + (- 1 ) m Q l2) h...i m (d il . . . d im 8){2-X) (12) 

where the multipoles are separated by X (see figure [ij. As our field theory is linear, 
the potential is simply the sum of the two multipole potentials 

cf>(x) = (-iyQW il ...i i (9 il ...a i ^)(x) + (-irQ^ 1 ... im (^ 1 ...a im ^)(f-X) . (13) 

Here (f(x) = — e~ m \ x \ j (47r|a?|) is the Green's function of the massive scalar field theory. 

The interaction energy is the difference of the total energy and the energies of the 
individual multipole fields, i.e. E int = E{<j>, p)-E(4> ( - 1 \ p^)-E{(jP\ p^>). Generally, 
the energy of point charges and all their derivatives is infinite, however, the difference 
-Bint is finite, as can be seen from the following equation which is easily derived from 



the original expression (eq. (10)) 



Efap) =E& 1 \pW) + E{<tP\p {i) ) 

3 x((9^( 1 ))(a^( 2 )) + TOV 1 )^ 2 ) + p( 1 )^ 2 ) + p( 2 )0( 1 )) . (i4) 

The last integral is therefore the interaction energy. Using the field equations, this 
can be reduced to 

£ in t,Yuk= fdW 1] <t> {2) ■ (15) 



We can as well write i?i n t,Yuk = / dPxp^tj)^ by exchanging the roles of the two fields. 
Using our expressions for p^ and 4>^ l \ we get 

£i„t,Yuk - (-l) l Q {2} i 1 ...i m Q {1) j 1 ...j l {d it . ■■d im dj 1 . ..dj,ip)(X) . (16) 

This is a surprisingly simple and symmetric result. In order to find the interaction 
energy between any two multipoles, we just need to calculate the appropriate deriva- 
tive of the potential (Green's) function. Note that compared to the interaction energy 
of the Coulomb potential 

iWouiomb = (-i) l+1 Q (2) i 1 ...i m QV jl ... jl (d il . ..d im o h . ..ay f --4=r j (it) 



X2 
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which is derived in eq. (8)], equation (16) here differs by a sign and the different 
choice of potential. The opposite sign is due to the Skyrme theory being a theory of 
scalar interactions, whereas electromagnetic interactions are vector interactions. 



2.1 Yukawa theory for asymptotic interactions 

In this section, we are going to show that to find asymptotic interactions of Skyrmions, 
it is sufficient to understand Yukawa multipole theory and know the falloff of the 
Skyrmions involved to get the complete picture. 

We compute the interaction energy using the product ansatz: given two Skyrmions 
defined by their Skyrme fields and U^ 2 \ we approximate their combined fields 
by U = U^U {2 \ The interaction energy is then defined by the equation 

E(U) = E(UW) + E(U&) + E int . (18) 

It is then useful to write the SU(2) fields in terms of their logarithms, i.e. 

U {1) {x) =exp (iF^^x)^ (19) 
U {2 \x) = exp(iF^ j {x)r j ) . (20) 

For convenience, we will assume our Skyrmions are centered at (0,0,0) for f/W 
and X = (X, 0, 0) for , where X is large. 

Let us write down the expression for the total energy E(U): We will immediately 
split up the energy into two integrals over the domains Q^ 1 ' — {x e ]R 3 |a;i < X/2} 
and ={i£ K 3 |a;i > X/2} and linearise = 1+iFj (x*)Tj where appropriate. 
Using integration by parts, this yields 

E{U) ^ E{U {1) ) + E(U {2) ) (21) 



12 
1 

12^2 
1 

12^2 



l — J d 3 xTr Uf^hth U + \[L^ 3 , [L^-, L< a >,]]) + m 2 U^ \ 

f dx 2 dx 3 Tr ( iFW„n ( L«i + -\L^ 3 , [L^^L^h]}) 

{X/2}xR 2 V V 4 / 

We can eliminate all the volume integrals using the field equations ([2]), and what 
remains is 

E int = / da; 2d^3 (FW&FWi - FM&FWi) (22) 

to linear order (this equation was first mentioned in the context of the Skyrme model 
in [10]). Doing the same computation in the massive scalar field theory defined 
by equation (111, i.e. computing the interaction energy of two charge distributions 



centered at (0,0,0) and (X, 0, 0), again using the field equations to reduce it to a 
planar integral, yields 

£int,Yuk = / dz a dx 3 (> 2 >cV (1) - (1) 9i^ (2) ) (23) 

J{X/2}xR 2 K ' 

with 0w generated by the charge distribution p( n \ 

This means to linear order (which is correct asymptotically, i.e. for X —¥ oo if the 
fields fall off quickly) the interactions are described by the massive scalar field theory. 
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Thus, given a 2 l - and a 2 m -pole in the Skyrme theory, their interaction for general 
separation X is 

Emt = g^(-l)'Q i ' (2) il ...i m Q i ' (1) J1 ... J1 (3 il ...d im d h ...d^)(X) (24) 



where the different prcfactor of the integrals in equations ( 22 ) and ( 23 ) has been taken 
into account, as well as that there are three scalar fields Fi in the Skyrme theory. 



3 The interactions of two B = 4 Skyrmions 

3.1 The B = 4 solution 

The Skyrmion solution with baryon number 4 (figure [3]) , first given in , plays an 
outstanding role in the Skyrme model. When quantised, it serves as the model for 
nuclei with four nucleons, notably the 4 He nucleus or a-particle, but also 4 H and 4 Li, 
which are however only seen as resonances and not as stable nuclei [T2"] . 

Interestingly, heavier nuclei seem to be composed mostly of interacting clusters 
of a-particles, a fact known to experimentalists and understood to a certain extent 
in the Skyrme model |13| . Nuclei with a baryon number divisible by four tend to be 
particularly well described by an arrangement of a-particles without any remaining 
single nucleons, and this is where the Skyrme model has been the most successful to 
date. This is why we are interested in knowing how the B = 4 solution interacts with 
copies of itself and other solutions of the Skyrme model. 

The B — 4 solution looks like a cube and has cubic symmetry. We will often 
be referring to the standard orientation of this solution in this paper. This is the 
orientation arising when one constructs the solution starting from the rational map 

R(z) = z4 + 2 ^ 2 + 1 (25) 
W z 4 - 2VSiz 2 + 1 

using the rational map formalism described in |14) . In this formalism, the rational 
map encodes the angular dependence of the Skyrme field, and an additional function 
f(r) gives the radial dependence. 

3.2 The symmetry group 

Let us describe the cubic symmetry of the B — 4 solution. The symmetry acts as the 
Oh subgroup of 0(3) on the spatial coordinates 12J, which is also known as the T± u 
representation of O^Q ^ does not act trivially on the Skyrme field. The action is 
easiest to describe as a group representation on the pion Fx, Fi and F3 fields. In the 
standard orientation the Oh symmetry can then be deduced from the rational map 



(eq. (25)) as acting by the E g representation on the (Fx,F%) field components and 
by A 2u on F 3 . 

There is a colour scheme corresponding particularly well to this symmetry, using 
Rungc's colour sphere to indicate the direction of the pion fields -7?/ j 7rj . The equator, 
the unit circle of the (tti,^) plane, is coloured using the hue colour attribute: This 
construction is also known as the colour circle. The tt 3 direction is indicated by light- 
ness, making the north pole (7r 3 = +1) white and the south pole (tt 3 = —1) black - 
see figure |4j In this colour scheme, the corners of the B = 4 cube in standard orienta- 
tion are black and white, whereas the faces are associated with the three fundamental 
colours green, blue and red. 



We are using the notation from 1151 for irreducible representations of O^. One-dimensional 
representations are labeled A% } A%, two-dimensional E and three-dimensional Ti, T2. Additionally, 
the index g (gerade=even) marks an even representation (where the inversion of Oh is mapped to the 
identity) and u (ungerade=odd) an odd representation (where it is mapped to the negative identity). 
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Figure 3: The B° = 0.1 isobaryon surface of the B = 4 solution in its standard orien- 
tation. The colours indicate the direction of the pion fields (see figure [2]). Opposite 
faces of the B = 4 solution have the same colour. 

3.3 B = 4 quadrupole and octupole tensors 

We want to apply the Yukawa multipole formalism to find the interactions of two 
B = 4 Skyrmions, so we need to identify the dominant pion multipoles of each of 
these. To this end, we are going to use the high symmetry of the Skyrmion to reduce 
the task as much as possible. 

We have three multipole expansions, one for each F^ field, whose multipole tensors 
we denote by Q\ i . We know how a multipole transforms under 0(3) as a tensor 
by applying the 0(3) matrix to each ij index. Oh, as a subgroup of 0(3), will act 
accordingly. However, there is another way to do the same transformation: doing the 
corresponding isospace transformation. Isospin transformations act on the i index. 
This gives us strong compatibility requirements, which we are going to explore next. 

As noted earlier, no static solution of the Skyrme field equation without pion mass 
can have a monopole. This is not true in the massive pion case. A monopole, being 
a scalar, will transform according to the A\ g representation of Oh- However [Fi^Fi) 
transforms according to E g , and F$ according to Ai u , so there is no monopole, even 
in the massive pion case. 

The next multipole to check is the dipole, which transforms according to the l u 



Figure 4: The colour sphere, used to colour the pion fields, in two different orientations 
(the second, top view shows the colour circle in the (^2,^2) -plane). An isobaryon 
surface of the B — 1 Skyrmion with coloured pion fields will have the same appearance 
as the colour sphere. 
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representation^] of 0(3), which corresponds to the T\ u vector representation of Oh- 
But under our isospace action, neither the .F^-fields nor Fg include any T\ u part, 
so the dipole can only be zero. 

The quadrupole, in our definition, transforms according to the 2 g © O g represen- 
tatiorj^jof 0(3), which corresponds to E g © T 2 g © Ai g of Oh- But we know Oh acts as 
on F3, so this field does not have a quadrupole moment. However, the (F±,F2) 
doublet transforms as E g , and E g is in the quadrupole representation - so it is possible 
to have a quadrupole in these fields. 

An explicit computation gives the following form for the quadrupole tensors: 



Q 1 = 4ttC 2 



Q z = 4vrC 2 -V3 . (26) 



This basis of traceless quadrupoles obviously depends on the orientation and is given in 
the standard orientation we are using. The single constant C 2 is left to be determined 
by numerics. 

The octupole transforms according to 3 U © 3 • l u in 0(3), which becomes Am © 
4Ti u © T 2u in Oh- This representation has no E g part, so the F\ and F2 octupoles 
are forced to zero. There is now, however, an A2 U part, so we can have an octupole 
in F 3 . 

This octupole tensor can be found to be 

Qw = 4 ^hik\ (27) 

where we write \£ijk\ for a tensor that is nonzero in the same components as the 
totally antisymmetric tensor, only that it takes the value 1 in all of these. 

3.4 The three-colour decomposition 

Looking at the B = 4 solution in figure [3] we see a very high degree of symmetry, re- 
flected also in the large symmetry group Oh- Despite this, at first sight the quadrupole 



tensors we found in equations ( 26 ) do not look very symmetrical in the axes. 

The reason for the apparent asymmetry of these tensors becomes clear when we 
look at the colour circle (figure [2]): In it, the three colours appearing on the faces of 
the B = 4 cube occur separated by angles of 2-7r/3 = 120°; however we are pressing 
them into the form of two orthogonal axes along F\ and F 2 . Would it be possible to 
rewrite the interaction in the form of three quadrupoles? The answer is yes, and the 
result is very neat and useful for us to get closer to a quadrupole interaction that is 
easy to work with. 



From equation ( 24 1 , we infer that the interaction of the two pairs of quadrupoles 
is given by 

El= 2 = Q^ijM^Q^u + Q 2 ' {1) tJ M ljkl Q 2 ^ kl (28) 

where we introduce the quadrupole interaction tensor = ^2(didjdkdiip)(X). 

Now we consider splitting Q 1 and Q 2 into three quadrupoles, one for each colour. 
Let us give the "green", "blue" and "red" quadrupole tensors numbers 1, 2 and 3 



3 We use £ g / u for representations of 0(3), where I is the angular momentum associated to the 
representation and the g/u index denotes even and odd, the same as used in Oh representations. 

4 It should be noted that this includes a trace part, whereas the common definition is for the 
quadrupole tensor in massless theories (e.g. electrodynamics) to be traceless. This however cannot 
be guaranteed in a massive field theory, and therefore we cannot a priori exclude trace parts in this 
and higher order multipole tensors. 
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Figure 5: The quadrupoles P 1 , P 2 and P 3 can be visualised as associated with vectors 
orthogonal to the faces of the corresponding colour. 



respectively and label these quadrupoles P c . Due to the 120° angles between colours, 
our decomposition should work like this: 

lp3 = Ql 

-P 2 = -io 1 -— O 2 
K 2 W 2 ^ 

where if is a constant. Note that P 1 + P 2 + P 3 = 0. We find that using if 2 = | the 
interaction energy in eq. ( 28 1 becomes 

E l=2 = p l,(l).. M . jkip l,(2) ki + p2,(l).. M .. kl p2,{2) ki + pUD ijMi]klP U2) u ( 30 ) 

This leads to an intuitive picture of one interaction term for each colour (green, blue 
and red). 

Assume for a moment that the Skyrmion is unrotated, giving Q 1 and Q 2 the form 
in eq. ( 26 ) . The three colour quadrupole tensors are 

/c l 

P 1 



p 2 = \/~4ttC 2 10 2 1 (31) 




P A = 



The form of these tensors is so simple that it lends itself to another interpretation: 
Note that the faces coloured green point in the x\ and — x\ direction. Blue has 
the same association with the a^-axis and red with the a;3-axis. Now, let us, for 
each colour, imagine a unit colour vector Sf pointing in the direction of one face 
of that colour (see figure [5| . It does not matter which of the two opposite faces we 
choose, as the vector will only occur in squared form. Then the corresponding colour's 

2 



quadrupole is given by P/i- = \ /|47rC2(3i^i^ — Jj,-). In the standard orientation, 



y V z ±\ % ] 

xf = Sir, so we can write 



P£ = d -AnC 2 (35 i£ S j£ - %) (32) 
using the convention that underlined indices are not summed over. 
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Now, making a spatial rotation £ SO(3) transforms P^ n ' to 

P 3 in) = ^ C ^ R f i 35 ^ '«) = yJl**Ck(MMli$ - %) (33) 
where we have reintroduced the n — 1,2 label for the Skyrmions at sites (1) and (2). 

3.5 Interaction of unrotated B = 4 Skyrmions 

In order to show how to use our three-colour decomposition to get concrete formulae 
for the interaction, we consider the case where RS 1 ' = — I. Also we are only 
considering the standard isospace orientation because in most, if not all of the known 
solutions composed of several B = 4 blocks, they are in the same isospace orienta- 
tion. This could be because changing isospace orientation would mix quadrupole and 
octupole interactions, likely leading to a weaker total interaction energy. 

In order to express the interaction energy easily, we need the n-th order modified 
spherical Bessel functions of the second kind, commonly denoted by k n {x) 1 which 
satisfy the differential equation 

! 2 ^P + 2x%M _ {x i + „(„ + l)) kn (x) = . (34) 
ax z ax 

The function k n decays to zero at infinity and has a pole at x = 0; it can be expressed 
in terms of elementary functions, using k (x) = e~ x /x and 

For convenience we use the rescaled functions 

(_ m \n+i 

Qn(r) = k n (mr) . (36) 

Then, the n-th derivative of the Green's function (p(x) — —e~ m \ x \ / (4tt\x\) can be 
written aa^ 

(0i 1 . . . 9i n tp) (x) — Q n ( |x | )i^ 1 Xi 2 ... Xi n -\- r^j [<^iiZ2 ^13 • • • ' ' ' ] 

\X\ 

~t~ |^| 2 [^ili2^3M^5^«6 ' " ' ^in ~t~ " " " ] "1" • * ■ (37) 

where successive terms contain one more 5-tensor. The notation \Ti^...i n 4- • • • ] de- 
notes the sum of all the non- equivalent terms obtained by permutation of the indices, 

e.g. [8 ll i 2 8i 3 i i H ] = 5^6^ + 6. lrl3 Si 2ii +5^5^. The quadrupole interaction 

tensor can then be expressed as 

M m = -L f q^X^XiXjXkXt + M^l^XfeX, + •••] + ijS kl + ■■■}) . 

6tt 2 y \X\ \X\ 2 J 

(38) 

Contracting the indices of Mijki with the quadrupole tensors P c as in eq. ( 30 ) , we find 
that the quadrupole interaction energy of two unrotated B — 4 Skyrmions separated 
by X reduces to 

<i 2 w = ^ (? 4 (m^ . (39) 



c=l 



5 Note that, if m = 0, this can be simplified using q„(r) = (-l) n+1 (2n - l)!!r-( n+1 '/(47r). 
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Now we turn to the octupole interaction. Using the octupole interaction tensor 



Mijklmn — ~7T^ { <?6 ( \X | )X i XjX k X l X m X n + ^ L - [Si j X k Xi X m X n 

6tH \ \X\ 



Q4(\X\) ~ - q 3 (|A|) \ 

Oijd k iX m X n H J H [oy d fc ;d mn H J (40) 



\X\ 2 J |X| 3 
we obtain the octupole interaction energy 



3 ^ ytHI 3 2|X| V 4j / |A| 2 |A| 3 / 

(41) 

The complete interaction energy of the Skyrmions up to octupole order is then 
E int 3 (^) = E \nt 2 (X) + E(= 3 (X). Note that the constants C 2 and C 3 depend on m. 
In the m = (massless pion) case, it is possible to simplify the interaction energy to 



Ant 



12Cf 1 



35^X^ + 21 



m=0 I I \ c=l / 

+ { A&2±1 ^ 1 + 21 E - 17 ) ■ ( 42 ) 

3.6 Interaction of B = 4 Skyrmions with 90° relative rotation 

Another interesting case is where one of the B — 4 Skyrmions is rotated by 90° 
relative to the other one around one of the cube axes. If the separation is along the 
same axis, this will produce maximum attraction for the octupoles, however slightly 
less attraction for the quadrupoles than in the unrotated case. This suggests it will 
be the preferred orientation at close range, but slightly less so at larger range. 

We choose the axis of rotation to be the xi-axis. The sense of the 90° rotation 
does not matter due to the symmetry of the solution. Using the same procedure as 
for the unrotated case but with 



R {2) = 1 , (43) 




we find the following result for the quadrupole interaction at general separation X: 
Ifif 12 

ElZot(X) ( «4(|*|)((3*? - If + 2(3X 2 2 - 1)(3X| - 1) 




(36X 2 -12) . (44) 



Symmetry tells us the octup ole part will change its sign compared to the unrotated 



case, « ot = ~E l = t 3 (eq. «4lb) and again e£*(X) = E l = t % t (X) + Ei=% t (X) 



As previously, a formula in terms of familiar functions can be obtained in the massless 
pion case: 



F i<3 

int,rot 



m=0 



— -i- (-35X? - 70X|X| + 20A 2 + 5 

TT \X\* V 23 

I- ( WlX\xlxl + 21 x c - 17 



15C 2 



\X\-r V 



(45) 



c=l 
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4 Numerical computation of the constants C2 and 

We have not yet determined the asymptotic coefficients Ci and C3, as they are not 
accessible to any but the crudest analytic approximations. To find C2 and C3, we 
need to solve a three-dimensional PDE, whereas the hedgehog asymptotic C\ could 



be found using ODE methods (see section 1.2 1 



4.1 Method 

A major problem with the computation of asymptotic values is the finite box size: 
Computing a solution of the Skyrme model in a finite box with vacuum U = 1 bound- 
ary conditions cuts off the asymptotics at the limits of the box. Worse, it has a large 
and hard-to-estimate influence even inside the box. Therefore, we choose a different 
approach: as the form of the asymptotic decay is known, we can impose this as a 
boundary condition, choosing some trial value for the constants C2,C3. Then the 
energy inside the box inside (C2, C3) can be computed using full three-dimensional 
numerics (by nonlinear PDE techniques), whereas the small energy outside the box 
£'outside(C2) C 3 ) can be estimated using the leading asymptotics by means of an inte- 
gration. 

The constants are then varied and the minimum of the total energy 
E = -Einside(C2, C3) + -Eoutside (C2 , C3) is determined. This will give their numerical 
values. 

To compute the energy inside the box of size L, we discretize the Skyrme en- 
ergy density using central 6th order finite differences on a cubic lattice with lattice 
spacing h = L/N, where N + 1 is the number of grid points in each direction. This 
gives the energy as a function of a discrete field. A nonlinear conjugate gradient 
methooj^] is then employed to find the minimum of this energy. The main program 
was implemented in the C programming language with thread-level parallelisation for 
computing gradients. 

To increase consistency, -Eoutsidc is determined using the same discretisation 
scheme as inside. The -Eoutsidc integral is computed in a box of total side length 16 
in the m = 1.0 case, 24 in the m — case considered below. These side lengths are 
chosen so that the neglected energy is lower than the error bound for the 
extrapolated total energy. We are first treating the m = 1.0 case interesting for 
nuclear physics. 



4.2 Numerical results 

It is expected that both the lattice spacing, h, and the box size L — N ■ h have an 
influence on the accuracy of the computations. Given infinite numerical precision, 
the limit of h — > and L — > oo should lead to convergence. Numerical precision and 
computation time put constraints on N however, which limits both L and h. In order 
to have a better control over the influence of h and L on the constants Ci and C3, we 
vary h and L determining the constants C2 and C3 in each case. 

For given h and L, E is computed for a grid of values of C2 and C 3 , leading to 
a graph as displayed in figure [6j A quadratic is fitted to these results to minimise E 
and arrive at the best estimate for C2 and C3 at this h and L. 

The different values are shown in tables Q] and [2] 

It can be seen from these tables that the influence of the lattice spacing h on the 
resultant parameters is rather small, whereas the box size L has a significant influence. 
To extract the final values, we take the h — 0.08 results and extrapolate the values 

^Precisely the "Nonlinear Conjugate Gradients with Secant and Polak-Ribiere" method described 
in |16| on page 53, without preconditioning. 
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Figure 6: The energy E for a Skyrme field inside a finite box can be computed for a 
range of values of C2 and C3. Taking the value of C2 and C3 where E is minimal leads 
to the static solution with the correct asymptotics. The crosses represent computed 
values, the mesh is a quadratic fit. This graph is for h — 0.08, L = 7.20. 



L 


1/L 


h = 0.2 


h = 0.16 


h = 0.1 


h = 0.08 


4.8 


0.208 


-1.783 


-1.779 


-1.766 


-1.761 


5.6 


0.179 


-1.841 




-1.838 


-1.836 


6.4 


0.156 


-1.870 


-1.873 


-1.871 


-1.870 


7.2 


0.139 


-1.888 




-1.891 


-1.890 



Table 1: Numerical value of C2 (quadrupole constant) depending on L and h at 
m = 1.0 



L 


1/L 


h = 0.2 


h = 0.16 


h = 0.1 


h = 0.08 


4.8 


0.208 


5.391 


5.389 


5.361 


5.350 


5.6 


0.179 


5.495 




5.501 


5.498 


6.4 


0.156 


5.531 


5.546 


5.546 


5.545 


7.2 


0.139 


5.548 




5.566 


5.566 



Table 2: Numerical value of C3 (octupole constant) depending on L and h at m = 1.0 



L 


1/L 


h = 0.2 


h = 0.16 


h = 0.1 


h = 0.08 


4.8 


0.208 


5.1794 


5.1808 


5.1815 


5.1816 


5.6 


0.179 


5.1786 




5.1802 


5.1802 


6.4 


0.156 


5.1784 


5.1796 


5.1800 


5.1800 


7.2 


0.139 


5.1784 




5.1799 


5.1800 



Table 3: Numerical value of E (the total energy of the B — 4 Skyrmion) depending 
on L and h at m = 1.0 
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Figure 7: C2 and C3 for different values of 1/L. The value at 1/L — > is determined 
by fitting a A + £? • {1/L) C function; A is then the extrapolated value. 



to 1/L — > using anA + B- (1/L) C ansatz (see figureM). Fitting these, we find 

C r 2 | m =i.o= -1.921 ±0.004 

C 3 |m=i.o = 5.581 ±0.003 . (46) 

Additionally, this technique gives us the possibility to get a value of the energy of the 
B = 4 Skyrmion at infinite box size. The same extrapolation technique applied to 
table [3] gives the value 

E\ m =x. = 5.1799 ± 0.0002 . (47) 

The errors given are taken from the fit errors of the extrapolation of 1/L — > 0, except 
for the energy, where the error from computing at finite h is dominant and is estimated 
using the results in table [3] 

The B = 4 energy at m = 1 has been computed before as E « 5.23 in [5] and 
E = 5.27 in |13) . While the deviation from our value seems significant, only a precision 
of about 0.5% was claimed there. 

The m = case is numerically somewhat more challenging, as the fields decay more 
slowly and the box has to be larger in order to obtain the same accuracy. From our 
experience in the massive case, we compute at different box sizes L for extrapolation, 
leaving h = 0.1 constant. Only for L = 8.0, we also compute at h = 0.16 in order 
to be able to estimate errors from the lattice spacing h (see tables |4]j6]). By fitting 
A + B ■ (1/L) C to the results in the h = 0.1 columns, we find the values 

C 2 \ m =o = -3.077 ±0.010 

C 3 | m=0 = 12.706 ± 0.020 (48) 
E\ m=0 = 4.4798 ± 0.0002 . 

The energy agrees with the value given in [U p. 377]. In this (m = 0) case, the fit 
error of the 1/L — > extrapolation is much smaller for both C2 and E, so that the 
error estimates above are the approximate variations between the h = 0.16 and the 
h = 0.1 values. For C3, both errors have a similar magnitude, so twice the fit error is 
given as the error estimate. 
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5 Application to B = 8, B = 12 and the Skyrme 
crystal 

5.1 The B = 8 Skyrmion 

A comparison of the asymptotic interaction energies of two B — 4 Skyrmions in the 
unrotated and 90° rotated case is shown in figure [8] This is for the separation vector 
X = (X, 0, 0). The asymptotic interaction energy in the unrotated case is slightly 
lower for large \X\. However, at short distances, the rotated configuration has lower 
energy, due to the large contribution from the negative octupole interaction energy. 
The latter going to — oo at the origin (zero separation) is clearly an artifact of the 
asymptotic approximation and the Skyrmions will repel at close distances. 



L 


l/L 


h = 0.16 


h = 0.1 


5.6 


0.179 




-2.506 


6.4 


0.156 




-2.763 


7.2 


0.139 




-2.892 


8.0 


0.125 


-2.968 


-2.961 



Table 4: Numerical C2 at m = 



L 


l/L 


h = 0.16 


h = 0.1 


5.6 


0.179 




10.960 


6.4 


0.156 




11.952 


7.2 


0.139 




12.352 


8.0 


0.125 


12.530 


12.520 



Table 5: Numerical C3 at m = 



L 


l/L 


h = 0.16 


h = 0.1 


5.6 


0.179 




4.4904 


6.4 


0.156 




4.4830 


7.2 


0.139 




4.4809 


8.0 


0.125 


4.4801 


4.4802 




Table 6: Numerical E at m = 



with 90 c 



unrotated 
relative rotation 



2.5 



3.5 



\X\ 



Figure 8: Interaction energies E^i for two B = 4 Skyrmions separated along the 
Xi-axis (m = 1.0); red/solid in the same orientation, green/dashed with 90° relative 
rotation around the xi-axis. While the asymptotic interaction energy in the rotated 
case goes to —00 as \X\ — > 0, there will actually be nonlinear corrections leading to a 
positive interaction energy at short distances. 
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For m = 1.0, the numerically determined B = 8 Skyrmion consists of two B = 4 
cubes, rotated by 90° around the axis of separation. In [17l p. 95], the distance 
between the cubes has been determined using the moment of inertia and found to 
be \X\ sa 3.0. At this separation, our calculations show clearly that it is energeti- 
cally favourable to rotate one Skyrmion by 90°, rather than leave it unrotatcd. The 
full interaction energy can be found using the energy of the B — 8 Skyrmion and 
subtracting two times the energy of the B = 4 Skyrmion (all energies taken from 
[13] ) and ii^j-Eint = —0.30. Inserting X = (3.0,0,0) into the asymptotic interaction 
energy, we predict a value of £^ 3 rot (A) = —0.42. This is surprisingly close to the 
full interaction energy given that the separation is relatively small. 



5.2 The linear B = 12 configuration 

One of the candidates for a minimal energy solution for B = 12 is a row of three B = 4 
Skyrmions, where the outer ones are in the same orientation while the middle one is 
rotated by 90° around their axis. Placing the three B = 4 cubes at —X = (— d, 0, 0), 
(0,0,0) and X = (d, 0,0), their asymptotic interaction energy is given by 

<-t 3 n„i2 = 2^ n l 3 rot(^) + ^ t 3 (2X) (49) 

where the interaction energies are as defined in sections |3.5| and |3.6[ Due to the 
exponential decay, the last term has a negligible influence, leading us to predict that 
d = 3.0 as in the B — 8 case (the same assumption was made in [T71 p. 106]). Here the 
predicted interaction energy is E il]tUnl2 = —0.85. Unfortunately, a numerical energy 
of the linear B = 12 solution of sufficient accuracy does not seem to be available at 
this point, so we cannot compare this prediction to any value in the literature. 



5.3 The Skyrme crystal 

The Skyrme crystal can be thought of as an infinite simple cubic (sc) lattice of B = 4 
Skyrmions, all in the same orientation (compare [18 ). Although we have found that 
at close range, a relative rotation of 90° will give a stronger face-to-face attraction, 
it is not possible to construct a crystal with all cubes in face contact rotated by 90° 
relatively around their separation vector. 

If we let a be the smallest distance between two B = 4 Skyrmions, then for any 
a we can compute the asymptotic interaction energy of an infinite lattice of these. 
Summing the nearest-neighbour contributions from the B — 4 Skyrmions adjacent 
through the 6 faces, the 12 edges and the 8 corners to a given B = 4 Skyrmion, and 
dividing by 2, this asymptotic interaction is 



pi<3 

int, crystal 



7 (128V3 - 3888 + 243^) Cf 1 



54tt 

m=0 



5 (4096^3 - 139968 + 28431 v^) Cf 1 



3888tt 



(50) 



per B = 4 Skyrmion for m = 0. In the massive pion case, the interaction energy can 
also be determined, and is illustrated below, but results in a far more complicated 
formula. 

The true energy E per unit cell in the crystal is rather easy to find numerically, 
as it just corresponds to finding a solution with B = 4 whilst imposing periodic 
boundary conditions. The box size L is then the same as the distance a between two 
neighbouring B = 4 Skyrmions. This has been done for the massive pion (figure [9]) 



and masslcss (figure 10 ) case. Additionally the predicted value from the asymptotic 



interaction is displayed. 

7 This number is the difference of two numerical energies of which the precision is not known and 
could therefore have a significant error. 
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Figure 9: The energy E per B ~ 4 Skyrmion for the Skyrme crystal at m = 1.0 
as function of the lattice constant a. The green/dashed line is the prediction from 
asymptotic quadrupole and octupole interactions. The relative error of the asymptotic 
approximation is also displayed. Note that for a > 7.5 the relative error is mainly 
due to the error in the total energy, suggesting the asymptotic interaction energy is 
actually much more precise. 



It can be seen that for m = 1.0, the asymptotic interaction energy predicts the 
crystal energy very well for lattice constants greater then about a = 4. The error 
however starts to grow again for large a. This is because it is numerically impossible 
to directly determine -Bint- To determine E lvi t, we have to compute the energy of the 
free B — 4 Skyrmion and subtract it from the energy in the unit cell of the crystal, 
which gives the determined E mt a large relative error. 

For the m = case, the prediction is not as good. This can be traced back to 
the lower precisions of the constants C2, C3 and E in this case (due to the higher 
computational complexity of determining them) and non- negligible contributions from 
more distant Skyrmions in the interactions. Also, due to the slower decay, it is 
expected that this interaction is only valid at much larger range than in the massive 
pion case. 

While in both cases, the octupole and the quadrupole parts of the interactions fea- 
ture different signs and the predicted asymptotic interaction therefore has a minimum 
(compare figure [8]), this is unfortunately clearly not in a range where the interactions 
can be trusted at all. Therefore, the asymptotic interactions could not be used to 
predict the Skyrme crystal lattice constant on their own: For m = 0, the lattice con- 
stant we find is a ~ 4.71 with an energy of E sa 4.15 per unit cell. The values found 
previously are a ~ 4.7 18j^]with an energy of E w 4.15 per B = 4 unit cell. In the 
m = 1.0 case, the Skyrmion crystal was not computed before. From the numerical 

8 Noticc that due to their definition of the Skyrme Lagrangian, length units are l/(eF 7v ) instead of 
2/(eF 7T ) used here. However, their L is the distance between two half-Skyrmions, and not two B = 4 
configurations as is our a here. So coincidentally the values of a here and L in |18| are numerically 
the same. 



18 



o 



s 



4.6 



4.5 



4.4 



4.3 



4.2 



4.1 



Numerically determined crystal energy + 
Relative error of asymptotic E- x - t x 

X 


h 

, t- +■ 

/+ 


1 1 + Jf-. 

x" 


_+ /+ 

/ + 
/ + 

X f 




/ + 




+ / + 




/ + 

+ / 

x 1 + 


X 


+ 1 + 




+ x :+Vx 
+ fx x 

++ x \ 


X 

X 


% XXXXX ^xxxxxxx v .., x X 





1.8 
1.6 
1.4 

1.2 fcf 



I 



1 

0.{ 
0.6 ^ 
0.4 
0.2 




10 



15 



20 



Figure 10: The energy E per B = 4 Skyrmion for the Skyrme crystal at m = 
as function of the lattice constant a. The green/dashed line is the prediction from 
asymptotic quadrupole and octupole interactions. The relative error for a > 15 is, 
similar to the m — 1.0 case, dominated by the error in determining E. For m = 0, 
we expect a larger relative error also due to contributions from more distant lattice 
sites as well as the interaction generally decaying more slowly. 



values presented in figure [9] here, we find 

a|m=i.o ~ 3.35 

E\ m=1 . w 4.95 per B — A Skyrmion (51) 

using a quadratic interpolation through the three values closest to the minimum. This 
can be compared to the energy value found in [TH] at m« 0.32, which i^E s» 4.3 at 
a 4.4. 

In figure |TT| two isobaryon surfaces of the m = 1.0 Skyrme crystal are displayed. 
Computing at the lattice spacing h — 0.1 restricts us to box sizes of a = (N + l)h, so 
we chose a — 3.30 which is very close to the minimum. 

For the m — Skyrme crystal, a translation by a/2 in one of the axis directions 
leaves the baryon and energy density invariant. Our isobaryon densities show that 
this is not the case for the m = 1.0 Skyrmion crystal, which is due to the asymmetry 
between U — 1 and U = — 1 created by the mass term (see eq. 0). In our config- 
uration, U = — 1 at the centre of the cube and U = 1 at the boundaries of the unit 
cell. 



6 Conclusion 

In this paper, we have shown in detail how to find the asymptotic interactions of the 
B = 4 Skyrmions, and laid out a general framework which can easily be applied to 
other solutions of the Skyrme model or soliton interactions in general, with massless 
or massive pions. 

9 Their L is given in the same units as used here, but a = 2L. 
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The B = 4 solution is very important in the Skyrme model, as it appears as a 
building block of many higher baryon number solutions. Working out its interactions 
can lead to a better understanding of how these B = 4 blocks arrange themselves in 
a more complex solution. With this work, the interactions have now been understood 
both qualitatively, by putting them into a nice computational framework using the 
three-colour scheme, and quantitatively, by determining the constants C2 and C3. 
This has been applied to the B — 8 and B = 12 Skyrmions to estimate their inter- 
action energies. A range of validity for the asymptotic interaction energy has been 
established using the Skyrme crystal. Together with C\, the hedgehog dipole con- 
stant, C2 and C3 encode our knowledge about the interactions of the most important 
solutions of the Skyrme model in the context of nuclear physics. 

Further research can now build on this framework to see how much the asymptotic 
interactions can serve to predict the form of solutions of the Skyrme model. It has 
been shown before that this approach works well in the case of four B = 1 hedgehogs 
making up a B = 4 solution, where the asymptotic interaction predicts how they 
should be oriented [B] . If the same could be done for solutions built from B = 4 blocks, 
or a mixture of B = 4 blocks and hedgehogs, this would constitute a tremendous help 
in the search for solutions. 

While the value of C\ on its own would have limited relevance to a qualitative 
structural understanding, as that is determined by the angular dependence of the 
interaction, the case of B = 4 is different (because there are two constants, so their 
relative strengths are important). In the case of systems consisting of B — 1 and 
B = 4 blocks, all three constants will come into play - and they could be interpreted 
as the structure constants of the Skyrme model, determining the composition rules 
of larger Skyrmions. 
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Figure 11: Two isobaryon (left: B° — 0.1, right: B° = 0.2) surfaces of the m = 1.0 
Skyrme crystal at the minimum (a = 3.3). It can be seen that a unit cell contains a 
cube that is more strongly connected internally than to the next unit cell. 
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